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Abstract 

In the context of the braneworld, a method to find consistent solu- 
tions to Einstein's field equations in the interior of a spherically sym- 
metric, static and non uniform stellar distribution with Weyl stresses is 
developed. This method, based in the fact that any braneworld stellar 
solution must have the general relativity solution as a limit, produces 
a constraint which reduces the degrees of freedom on the brane. Hence 
the non locality and non closure of the braneworld equations can be 
overcome. The constraint found is physically interpreted as a neces- 
sary condition to regain general relativity, and a particular solution 
for it is used to find an exact and physically acceptable analytical in- 
ternal solution to no- uniform stellar distributions on the brane. It is 
shown that such an exact solution is possible due to the fact that bulk 
corrections to pressure, density and a metric component are a null 
source of anisotropic effects on the brane. A conjecture is proposed 
about the possibility of finding physically relevant exact solutions to 
non-uniform stellar distributions on the brane. 
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1 Introduction 



The standard model and the general relativity represents the two great the- 
ories in fundamental physics. The success of general relativity is beyond any 
doubt, however due to its inconsistency with quantum mechanics, it is not 
possible to ensure that this theory keeps its original structure at high ener- 
gies. Indeed it is expected that such a structure will be modified and thus 
general relativity will emerge as the low energy limit of a more general theory 
not yet discovered. In spite of an enormous effort, it has not been possible to 
construct an unified theory which describes all fundamental interactions. To 
carry out this project, a consistent quantum theory to gravity must be cons- 
tructed. Fortunately there is a theory which seems to lead to the expected 
quantum gravity. 

The superstring/M-theory is considered one of the most promising candi- 
date theories of quantum gravity. It describes gravity as a high dimensional 
interaction which becomes effectively four dimensional at low enough ener- 
gies. This theory has inspired the construction of braneworld models, in 
which the standard model gauge fields are confined to our observable uni- 
verse (the brane), while gravity propagates in all spatial dimensions (the 
bulk) p. 

The implications of the braneworld theory on general relativity have been 
extensively investigated [2]-[ll] (see also a review paper |T5] and references 
therein), most of theme on cosmological scenarios [16]- (see also [27] and 
references therein). The studies on astrophysics consequences |28]-[36] are 
mostly limited to exterior solutions, even though it is well known that grav- 
itational collapse could produce very high energies, where the braneworld 
corrections to general relativity would become significant |37]-[l3]. Unfortu- 
nately there has not been reported any consistent internal stellar solution 0, 
therefore the braneworld consequences on physical processes in the interior 
remain unknown so far. 

As is well known, it is extremely complicated to find exact interior so- 

^In the pioneer work of Germani and Maartens [28] exact solutions for a uniform 
distribution were found, where it can be seen how much more difficult it would be to find 
a solution except for uniform stellar distributions. 
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lutions in general relativity [H]. Even for the simple case of a static per- 
fect fluid, only very few reported exact solutions are of physical interest 
[45j . Hence it is not surprising that finding exact solutions in the context of 
braneworld, where new terms from the bulk have to be considered, becomes a 
very complicated process, except when uniform distributions are considered. 

Consequently important issues, as the role played by the density gradients 
as a source of Weyl stresses, and its importance in different physical pro- 
cesses in the interior, have not been elucidated yet. An exact solution would 
be useful to carry out an analytic analysis. Unfortunately there has not been 
possible to construct any exact internal solution to non-uniform stellar distri- 
butions in the context of the braneworld. The reason is that the non locality 
and non closure of the braneworld equations, produced by the projection 
of the bulk Weyl tensor on the brane, lead to a very complicated system 
of equations which difficult the study of non-uniform distributions. Indeed, 
there is not any well established criteria about what restriction should be 
considered on braneworld equations to obtain a closed system To solve 
this issue, it is necessary a better understanding of the bulk geometry and 
how our 4D spacetime is embedded. 

As already mentioned, the study on density gradient effects in the interior 
remain unknown so far. For instance, the role played by density gradients 
during the gravitational collapse represents an open issue which has not been 
studied in the context of braneworld. In this respect, it is well known that the 
collapse of a homogeneous star leads to a a non-static exterior solution |38j - 
[H]. Moreover, Govender and Dadhich [IQ] proved that a collapsing sphere 
on the brane radiates, and they shown that the exterior for this radiative 
sphere can be described by a Vaidya metric that envelops the collapsing 
region. This represents a new and complete different aspect from the four- 
dimensional general relativity situation, and a remarkable issue which should 
be deeper investigated. In this context, it would be interesting to elucidate 
the effects produced by density gradients as a source of Weyl stresses, which 
eventually will be, together with the non- vanishing of the effective pressure 

^Recently some braneworld wormhole solutions were reported |46j . where an exact dust 
solution can be clearly seen. 

■^In [47j a closed system of equations for static spherically symmetric situation in the 
case oiP = was studied. 
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at the stellar surface, the responsible for the non-static nature of the exterior 
solution. 

In this paper we develop a method which helps to overcome the non 
locality and non closure of the braneworld equations when a non-uniform 
stellar distribution is considered. This approach is based in the fact that 
any braneworld stellar solution must have the general relativity solution as 
a limit. This statement is exploited in this paper to generate a constraint 
which helps to reduce the degrees of freedom for an non-uniform distribution 
having local and non-local bulk terms on the brane. A particular solution 
for this constraint is used to find an exact and physically acceptable internal 
solution to no-uniform stellar distributions on the brane. 

This paper is organized as follows. In Section 2 the Einstein's field equa- 
tions and matching conditions in the brane for a spherically symmetric dis- 
tribution is reminded. In Section 3 a toy solution is used to show how the 
general relativity limit could be easily lost when a solution is achieved on 
the brane. In Section 4 the source of the general relativity limit problem 
is clearly identified and a solution for this is given. This solution shows an 
explicit constraint emerging as a natural consequence of the general relativ- 
ity solution. In Section 5 the constraint identified in the previous section is 
used to find an exact analytical solution to a no-uniform distribution with 
bulk terms on the brane. In Section 6 an exact and physically acceptable 
solution is found. It is shown that this exact solution can be obtained due 
to the null consequence of bulk corrections to pressure, density and a metric 
component (p, p and v) to produce anisotropic effects on the brane. In the 
last section the conclusions are presented and a conjecture is proposed about 
the possibility of finding physically relevant exact solutions to non-uniform 
stellar distributions on the brane. 
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2 The field equations and matching condi- 
tions 



The Einstein field equations on the brane may be written as a modification 
of the standard field equations [lO] 



G 



(1) 



where A is the cosmological constant on the brane. The energy-momentum 
tensor has new terms carrying bulk effects onto the brane: 



T, 
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here a is the brane tension. The new terms 5*^1, and S^^, are the high-energy 
corrections and KK corrections respectively, and are given by 
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being U the dark radiation and V^u and the anisotropic stress and energy 
flux respectively. 



We consider a spherically symmetric static distribution, hence = Q 



and 



(5) 



where is a unit radial vector and h^^ = g^^ — u^Uu the projection tensor 
with 4- velocity u^. The line element is given in Schwarzschild-like coordinates 

by 

ds^ = e^dt^ - e^dr^ - (dO^ + sin ^Odcf) (6) 
where v and A are functions of r. 

The metric ([6]) has to satisfy ([1]). In our case with A = we have: 



1 + fl-^ 

1^2 I ^2 ly' 
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where /i = df/dr and fc^ = Svr. The general relativity is regained when 
— > and ffTOj) becomes a lineal combination of ©-Q. 

The Israel-Darmois matching conditions at the stellar surface S give 







where = /(r) —f{r) Using ffTTl) and the field equation ([T]) with 
A = we have 



which leads to 



(12) 



0. 



(13) 



This takes the final form 



1 
cr 



Pi 
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where //j = /(r) \r=R- The equation gives the general matching condi- 
tion for any static spherical star on braneworld [28] @. When we 
obtain the well known matching condition p/j = 0. In the particular case of 



^The general matching conditions on the brane for a spherically symmetric vacuum 
region embedded into a cosmological environment can be seen in j48|-[50j 
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the Schwarzschild exterior solution lA^ = = 0, the matching condition 
f fT^ becomes: 

i (I + + ^51.0^ (15) 

It is easily seen that the matching conditions do not have a unique solution 
on the brane. 



3 Losing general relativity 



In the case of a non-uniform static distribution with local and non-local bulk 
terms, we have an indefinite system of equations (|71)- (fT0l) whose solution must 
be determined with the add of additional information. However it is not clear 
what kind of restriction should be considered. As mentined earlier, this is an 
open problem which solution imply more information of the five dimensional 
geometry and a better understanding of how the brane is embedded in the 
bulk. Hence the study of any method which leads to a reduction of the 
degrees of freedom on the brane represents a subject of great interest. In 
order to construct a method, let us start from the simplest way to find a 
solution to the system fl7j)- f|T0l) 

From the field equations ([8]) and iQ we obtain 



and 



with 
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Now using (fT7|) in the field equation ([7]) we have 



A I „-A 



Vn + z/2/2 + 2z/i/r + 2/r2 



1^1 /2 + 2/r 



r2(z/i/2 + 2/r) 



-»7r 



(z/i/2 + 2/r) 



(20) 



which formal solution is 
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87r(p -3p (p^ + ^PP) 

a ^ 



rfr + c , (21) 



'-i (f + f) ^''^ 

Hence finding z/, p and p satisfying ffTOj) we would be able to find A, V and 
U by f l2Tl) . ( fT6l) and ( fTTl) respectively. Therefore finding a solution to the 
system (171)-( |T0|) . at least from the mathematical point of view, seems not 
very complicated. We will see that even from a pure mathematical point 
of view, to find a consistent solution is not easy at all. To address this we 
consider a simple solution to (|TOl) showed below 

p = A + 5r^/2. p^_^_|^A./2. ^u/2^^^k/2^ (23) 

where > 0, A, i? and C are constants. Although such a solution is not 
of physical interest, it will be useful in showing the problem which arises 
when the complete solution is found by this method. In this sense, (!23l) is 
considered to be just a toy solution. 

Using (123!) in (pT!) we should expect to find a solution given by 

e"^ = 1 - — /' r^pdr + -(Bulk effects), (24) 
r Jo cr 

and then taking the ^ — limit the general relativity would be regained. 
Unfortunately this does not happen, thus the solution found by the simple 
way showed here is not even mathematically consistent, that is the solution 



8 



evaluated at ^ = is not a solution to ©-(IHl) at ^ = 0. Hence the appar- 
ently simplest and straightforward way lead us to an inconsistent solution, 
so obviously there is something incorrect in the way we just used to obtain 
A(r), therefore it is necessary to carry out a more careful analysis. This will 
be addressed in the next section. 



4 Recovering the general relativity limit 



In order to find the source of the "general relativity limit problem" found in 
the previous section, let us write the field equation ([7]) as usual in general 
relativity 

e-^ = l- — TrV^, (25) 
r Jo 

where 

(2a, 

Thus we can see the high-energy effects of bulk gravity and nonlocal correc- 
tions from the bulk Weyl curvature. We can see that the "solution" (l25il 
depends itself on A(r) and Ai(r) through U (fT7|) - (fT9|) . hence it it represents 
an integral differential equation for A(r). Now we can realize why the formal 
solution (12T]) cannot written by the way showed in ( |24l) . Such solution has 
mixed up general relativity terms with bulk effects on the brane. In conse- 
quence finding how the general relativity limit problem originates becomes 
very complicated. To clarify this point let us rewrite the differential equation 

m by 



-Aie"^ e-^ 1 
\ 5 5- + 87rp 



-A 



+ 2— 



SttSp p (p + 3p) 

cr 



0, 



(27) 



here the left bracket has the standard general relativity terms and the right 
one the high energy terms with bulk Weyl curvature contribution to the 
differential equation of A(r). It is worth noticing that in the right bracket 
only high energy terms are manifestly bulk contributions, hence to keep bulk 
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Weyl curvature contributions under control when the formal solution fl2ip is 
achieved is not at all easy. 

We propose a solution to fl27p which can be written in a way that allows 
to see clearly the effects of bulk contribution on the brane 



8^ r 2 , -/ r 

1 r pdr + e — 

r Jo JO Mtt 



2 ~ r' 



6n 



H{p,p,iy) + — [P+ 5pp 



where / is given once again by fl22|) and 

H{p,p, u) = 8n3P 



dr, 
(28) 



,z/i 1, , z/? 2i/i 1 , 

+ -) + f^i^u + V + — + 
/ r z r r 



with 



/i = 1 / r^pdr. 

r Jo 

The reader may check that ( 128|) satisfies (1271) . 

It is easy to see by ©-(in]) that at = we have 



, (29) 



(30) 



2Gl + Gl 



,1/1 1, , z/? 2z/i 1 

f^ii^ + -)+ /^(^ii + T + — + - 
2 r 2 r r 



SrrSP, (31) 



so the function H{p^p,v) vanishes when cr^^ and thus it may be in- 
terpreted as a function which measures the anisotropic effects due to bulk 
consequences to p, p and z/. Since H{p, p, z/) vanishes when the 

solution fl28l) seems having the well known general relativity limit: 



r 2 . 

— r pdr. 
r Jo 



(32) 



However using our toy solution fl23l) we can see that H{p, p, v) does not van- 
ish. We should expect this to be due to the fact that this solution was built 
using only flTUl) . without any additional assumption. Hence we may learn 
that there is something else to be considered to find at least a mathemati- 
cally consistent solution. In order to close the system (17|l- (fT0|) . at least two 
conditions must be imposed. Thus we could use if(p, p, z/) = as a conve- 
nient constraint which lets us make sure that A(r) has the general relativity 
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solution as a limit. In any case every solution found must satisfy the condi- 
tion H{p, p,!/) when 0. In that sense this limit could be thought 
of as a constraint itself which might be useful to close the system. However 
we cannot forget that the H{p, p, v) contribution to A(r) is under integration, 
therefore the constraint to be considered is 

lim^-i^ / ^ H{p,p, iy)dr = 0, (33) 

and it must be satisfied as a necessary condition to regain general relativity. 
Since the constraint fl33l) is proportional to H{p,p, u), we may interpret this 
constraint as the continuous elimination of deformations produced on general 
relativity by anisotropic effects generated for bulk corrections to p, p and u. 

The method shown here leaves us with three unknowns p, p and u sa- 
tisfying (fTOl) . and eventually the constraint (133|1 . However there is something 
which we might be concerned about. We can see that Eq. (fTOl) does not have 
any manifestly bulk contribution. Thus it might be thought that the solution 
eventually found for p, p and u would be the same as the general relativity 
one. Hence there would not be any bulk contribution for these functions. 
Actually the fact that Eq. ( ITOi) does not have any manifestly bulk terms does 
not mean its solution does not either. Indeed the bulk contribution for p, p 
and u will appear by means of matching conditions, where the assumption 
of vanishing pressure will be dropped p9], |13] as will be shown in the next 
section. 



5 A solution 

We shall construct a simple mathematically consistent solution in order to 
clarify the method described in the previous section. First of all let us enforce 
the constraint 

H{p,p,u) = (34) 

on the brane to make sure we have a solution for the geometric function A(r) 
with the correct limit. Of course satisfies (l33l) and what it means is that 
eventual bulk corrections to p, p and z/ will not produce anisotropic effects 
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on the brane. Secondly, to evaluate (12^ and eventually (1251) . we need to 
consider a simple enough expression to z/(r). Accordingly we impose 

e'' = Ar\ (35) 
A simple but unique solution to ffTOj) and flM|) using fl35l) is given by 

p(r) = 5r-^/3. 8^p(^) = ^ - 87^35r-^/^ (36) 

where A and B are constants to be determined by matching conditions. 
Using ([35]) and ([Ml) in ([28]) we obtain 

. 2M/r\5/3 SttE /18 ./o 125 ./.A , , 

where 

M= / Anr'^pdr, (38) 

JO 

and where R is the radius of the distribution. We see that the solution thus 
generated is not regular at the origin, caused by the imposed function ( l35il 
and the particular solution (136!) . but not due to the method described here. 
Finally using (fT6l) and ( fTTj) we obtain 

4^ = ^ (^;:^ - Y?^) ; ^ = " Ir^) • ^^^^ 

The bulk contribution to p, p and v can be found by matching conditions. 
For instance let us consider the Schwarzschild exterior solution 

e^^ = e-^+ = 1 - W+ = P+ = 0. (40) 

r 

Considering the matching condition = at the stellar surface S we 

have 

AR^-l-'-^. (41) 

2A4 ^ 2M _ SttE /18 4/3 _ 1^ o-2/3\ .43^ 
i? i? a U3 7 / ^ ' 
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and using (fT^ we obtain 



B{a) 



^765 - 26 5 27ri?2 + VsT J^ 11475 + i286007ri?2 + I379047r2i?4 



i62407ri?2/3 



(43) 

thus the bulk contribution to p, p and v can be seen clearly through A{a) 
and B{a). Considering lineal terms in we have 



B(a) 



1 



+ 



1 



245 



67ri?2/3 a 159127r2i?8/3 



1 4 



M(a) = M 

^ ^ aR V13 2l7r 



yl(a) 



1 



2M 1 



1 



ai?2 Vl3 2l7r 



(44) 
(45) 
(46) 



Thus we obtain our solution at first order in 1/cr, where the bulk effects on 
p, p and z/ can be easily seen by 



P 



57r 



-3 



1 1 

+ 



245 



67ri?2/3 (T 159127r2i?8/3 
1 245 



-4/3 



+ 



67ri?2/3 a 159127r2i?8/3 

2M 18/3 1 
1 - — + 



ai?2 Vl3 2l7r 



4/3 



leaving 



7]-2/^2/3^8/3 
1 



4P 

cr 

w 

a ~ 487r2i?2/3^8/3 



52r2/3 153i?2/3 
1 1 



and 



2M/r\5/3 1 



13r2/3 17i?2/3 
8 / 3 1 



1 

a' 
1 

0" 



1 



ai?2/3^2/3 Vl3r2/3 2l7ri?2/3 



(47) 

(48) 
(49) 

(50) 
(51) 

(52) 



The solution P71) - fl52l) represents an exact analytical solution to the sys- 
tem fl71)- f|T0l) with the correct limit at low energies. The algorithm developed 
to obtain this solution runs as follows: 
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• Step 1: Impose the constraint H{p^p,v) = to make sure we have a 
solution for the geometric function A(r) with the correct hmit at low 
energies. 

• Step 2: Pick a simple enough expression to z/(r) to obtain an analytical 
expression to the Eq. ( I22l) . 

• Step 3: Use v{r) in both the conservation equation p' = — y(p + p) 
and the constraint H{p,p, u) = to find p{r) and p(r). 

• Step 4: Use p(r), p(r) and z/(r) in (!28l) to obtain A(r). 

• Step 5: Find V and U using (fT6|) and (fTTl). 

• Step 6: Drop out the condition of vanishing pressure at the surface to 
obtain the bulk effect on any constant C C{(^)- Then we are able to 
find the effects of bulk gravity on pressure and density. 

It is worth noticing that this algorithm is general and works even for a 
numerical analysis, where any arbitrary expression for z/(r) at step 2 can be 
chosen. 

So far the solution found is a simple mathematically consistent solution. 
It was constructed to clarify the method described in the previous section 
and summarized in the algorithm shown here. In order to obtain a physically 
acceptable solution, it is necessary to carry out a more careful analysis. This 
wiU be addressed in the next section. 

6 A physically acceptable solution 

The approach developed here ensures general relativity as a limit, but it does 
not say anything about finding a physically relevant solution. Of course it 
is not enough to pick out a simple expression to z/(r). Indeed choosing a 
convenient expression for z/(r), leading to a physically acceptable solution, 
represents a critical issue. In this paper it was found that the following 
expression for the metric component 

e" = A(l + Cr2)3 (53) 
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leads to a simple and physically acceptable solution for (ITU]) and given 
by 

^ 9C(1 -Cr2) 

Svrpfr) = ^ ^, (54) 

and 

3C(3 + Cr^) 
87rp(r) = ^, (55) 

where A and C are constants to be determined by matching conditions. 



The solution for the geometric function A(r) is obtained using (!53l) -(l54 
in ( |28l) . leading to 

(56) 



e-A{r) ^ 1 _ 2m(r) ^ 
r ' 

where the interior mass function m is given by 

1 g{r) 



m[r] = m[r] 



a 167r(l + Cr2)(2 + 5Cr2)ii/io 
with 



(57) 



. ^ nr^2 rr2(2 + 5Cr2)T^ (_9 + 3C7r2 + 2^^^) , , ^, 

-'0 (l+Cr^) 

and m(r) being the general relativity interior mass function, given by the 
standard form 

m(r) = AT^r pdr = j^^^^^, (59) 
hence the total general relativity mass is obtained 

where i? is the radius of the distribution. A numerical analysis on g{r) shows 
that g{r) > for all r < R. Hence by (IFH) we can see that the effective 
interior mass function m(r) is reduced by five dimensional effects. 

Using (1161) and (fT7|) the interior Weyl functions are written as 

, 2 



P(r) 



4r3(l + Cr2)^(2 + 5Cr2)T5 



2(l + Cr2) (l + 7Cr2 + 14CV^)^(r) 



+3 C (2 + 5 C ry " 18 - 1 1 1 C - 137 + 86 C + 40 

(61) 
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U{r)= ^ -U(i + Cr^)' (b + lCr^)g{ 

16r(l + Cr2)'(2 + 5Cr2)Ts V ^ ^ ^ ^ 



+ 3 C r (2 + 5 C r^) '° (-180 - 660 C - 509 + 62 + 55 

(62) 

It can be seen that A, U and V do not have exact expressions due to the 
presence of g(r). However this feature does not impede the construction of a 
physically relevant exact solution for the pressure and density, whose explicit 
brane world expressions can be seen using matching conditions, where the 
assumption of vanishing pressure will be dropped [29], |13]. As Schwarzschild 
is not the only possible static exterior solution, we have many scenarios to 
consider. Let us begin considering the simple Schwarzschild exterior solution 

e^^ = e-^^ = 1 - W+ = P+ = 0. (63) 

r 

The matching condition [ds^Js = at the stellar surface S yields 

A = (l-^)(l + Ci?V, (64) 



2A1 ^ 2M_1 g{C) 

R R (T87ri?(l + Ci?2)(2 + 5Ci?2)ii/io' ^ > 

where 

9{C) = g{r)\r=R = ~9C (l + Cr^f ^'^ ^ 

Using (|T5|l it is found that C must satisfy the condition 

1 9iC) (l + 7CR^) 
SnpiR) + ^^^^ = 0, 67 

^ 8 TT i?3 (1 + C (2 + 5 C i?2) To 

which can be written as 

9C{1-CR') 1 g{C) {1 + 7CR') 

^ ~ inr ~ ^- \^°) 



2{1 + CR^y f^8 7ri?3(l + Ci?2)2(2 + 5Ci?2) 

16 
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The equation fl68l) shows clearly that C cannot be the general relativity value 
Co, given by| 

Co = (69) 

which is found using the condition p{R) = in (l54l) . What the equation (168!) 
says is that the general relativity value of C, namely Cq, has been modified 
due to the buck effects by 

C = Co + 6ia). (70) 

In this sense S{a) represents the " bulk perturbation" of the general relativity 
value of C. In order to find the bulk contribution to p and p we need to find 
C, given by ( |70l) . satisfying ( |68l) . Thus we have at first order in 

The pressure can thus be found expanding p{C) around Cq 

p{Co + 6)=p{Co)+6^\c=Co, (72) 

which leads to 

9Co(l-Co r^) 9 1 - 3C,r^ 

5p((T) 

and by the same way the density is found to be 

3 Co (3 + Cor^) 3(3 -Cor^) 

167r(l + Cor2) ^167r(l + Cpr^)^ 

5p(a) 

where Sp{a) and 6 p{a) represent the bulk effects on p{r) and p(r) respec- 
tively. However at the surface and for any arbitrary R always we have 

PiR) = --^^i^) < 0- (75) 



numerical analysis shows that g{C) 7^ for a wide range of R. 
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Hence the Schwarzschild exterior solution is incompatible with the interior 
solution found here. Thus a different exterior solution must be considered. 

Using now the Reissner-Nordstrom-like solution given in [3] 

e'* = e-* = 1-^ + 4, (76) 

and considering the matching condition [ds'^]^ = at the stellar surface S, 
we have 

A(l + Ci?T = l-^ + Jj, (78) 

2A1 ^ 2M_1 giC) a 

R R a87ri?(l + Ci?2)(2 + 5Ci?2)ii/io ^ i?2' ^ ) 

and using (IT^ we obtain 



2 ^ 9Ci?^(l-Ci?^) 1 il + 7CR')giC) 

2(1 + Ci?2)2 a87ri?(l + Ci?2)2(2 + 5Ci?2)- ^ 

The constants Ai and g are given in terms of C through equations (1791) 
and ( l80l) respectively, leaving A and C satisfying ( 178|) . As in the previous 
case, C and A have a well definite general relativity values, which are given 
respectively by the condition p(i?) = in (15^ and evaluating (ITHIl at a^^ = 

Ao(l + Coi?^)3 = l-^. (81) 

JrC 

Comparing ( 178|) from ( IHTi) it is clear that the general relativity values of A 
and C have been modified by bulk effects by 

da) =Co + 6iay, A{a) = Ao + eia). (82) 
Using (El in ([78]) we obtain 

{Ao + e)[l + {Co + 6)Rr = 1 - ^ + J^- (§3) 
18 



Evaluating the expressions (17^ and 

the equation fl83l) leads to 



terms in a 



dX C = Cq + 5 and keeping lineal 



C{a) = Co + 



1 



a 167ri?(7) 
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(84) 



Hence C(cr) is determinate if e is kept as a free parameter which can be 
used in finding a physically acceptable model. Thus it is possible to see bulk 
consequences on p and p through (173|) and (17^ . The figure 1 shows the 
behaviour of the pressure in both the general relativity and braneworld case. 
It can be seen that the five dimensional gravity effects reduce the pressure 
deep inside the distribution. However the situation changes for the exterior 
layers, where the matching conditions lead to p 7^ at the surface. 



7 Conclusions 

In the context of the astrophysics braneworld, a method to find consistent so- 
lutions to Einstein's field equations in the interior of a spherically symmetric, 
static and non uniform stellar distribution with Weyl stresses was developed. 
During this process a very important feature of the indefinite system ([7])- 
(fTOj) was identified. When the simplest and straightforward way to solve 
field equations is used, the general relativity solution is lost. The source of 
this problem was clearly identified, and a general solution for it was given 
through the equation fl28|) . which gives the geometric function A(r) written 
in a way that allows to see clearly the bulk effects on the brane. 

Keeping under control non-manifest bulk contributions when the solu- 
tion for A(r) is achieved on the brane, the loss of the general relativity limit 
was avoided and then a new constraint arose as a natural consequence of 
this method. A physical interpretation of this constraint was given as the 
necessary condition to regain general relativity. It was shown that this con- 
straint can be used as a natural tool which helps to close the indefinite system 
fl7j)- f|T0l) and therefore to investigate bulk consequences on the brane. 

An exact analytical solution on the brane was constructed by a particular 
solution to the constraint (|33|1 . This particular solution represents a new 
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constraint defined on the brane, which was interpreted physically as the null 
consequence of bulk corrections to p, p and v to produce anisotropic effects 
on the brane. 

By prescribing the temporal metric component Qoq, an exact interior so- 
lution to the pressure and density was found. This represents a physically 
acceptable solution, namely, regular at the origin, pressure and density defi- 
nite positive, well definite mass and radius, monotonic decrease of the density 
and pressure with increasing radius, etc. It was shown that this solution is 
incompatible with the Schwarzschild's exterior solution. Using the Reissner- 
Nordstrom-like metric given in [1], the effects of five dimensional gravity on 
pressure and density were found through matching conditions. It was found 
that the bulk gravity effect reduces the pressure deep inside the distribution, 
but the situation changes for the exterior layers, as a direct consequence of 
matching conditions. 

As remarked earlier, it is well known that finding physically relevant exact 
stellar solutions in general relativity, even for a perfect fluid, is not at all easy. 
Due to local and non-local contributions from the bulk, the problem becomes 
even more complicated. The exact solution for the physical variables p{r) and 
p(r) found in this paper was possible as a direct consequence of the constraint 
It is easy to see through the field equations ([7]) -([9]) evaluated at cr~^ = 
(general relativity) that the non-local function H{p,p, u) can be written as 

H{p, p, u) = [2Gl + Gl) U-i=o - 87r3p, (85) 

which clearly correspond to an anisotropic term. In general relativity the 
function H{p, p, z/) vanishes as a consequence of the isotropy of the solution. 
However, in the braneworld case, the condition H{p,p,v) = in general is 
not satisfied anymore. There is not reason to believe that the modifications 
undergone by p, p and due to five dimensional effects on the brane, do not 
modify the isotropic condition H{p, p, u) = 0. Therefore in general we have 
H{p,p, u) ^ on the brane. However when the constraint H{p,p, z/) = is 
imposed to ensure the general relativity limit, all anisotropic effects on the 
brane whose source are bulk corrections on pressure, density and u, are elim- 
inated, leaving only high energy corrections for the geometric function A(r). 
Such constraint not only ensures general relativity as a limit, but also rep- 
resents an enormous simplification for the function A(r), which clearly helps 
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when exact solutions are investigated. Hence this constraint should be con- 
sidered in searching physically relevant exact solution on braneworld. Guided 
by this, the following conjecture is made: When a spherically symmetric, 
static and non-uniform stellar distribution in the context of the braneworld 
is considered, it is not possible to generate physically relevant exact analytical 
solutions on the brane, unless bulk corrections on pressure, density and the 
temporal metric component do not produce anisotropic consequences. This 
conjecture might help in the search of exact and physically relevant solutions 
when a non-uniform distribution is considered. Hence the role played by 
the inhomogeneousness as a source of Weyl stresses in the interior could be 
studied. This is currently been investigated. 
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